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ABSTRACT 


A time-dependent, two-dimensional thermodynamic model of an open lead 
in central arctic sea ice is presented. The model is generated by opening 
a lead of finite width and infinite extent in the equilibrium sea-ice cover. 
From this initial condition, the model is integrated numerically as a sea- 
ice cover is reestablished over the lead. The effects of various representa- 
tive advective parameterizations, and temperature and salinity profiles in 
the ocean's surface layer, on the heat flux through the lead and on the 
nature of ice formation are investigated as the lead closes by thermal 
processes alone. Continuity equations involving a horizontal advection 
term and a vertical diffusion term govern heat and salt transport in the 
water. Cooling-induced convective overturn as a mechanism for vertical 
neat anda salt transport in the water column is treated through an artifice 


of the vertical diffusion term. 
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I. INTRODUCTION 


A. POLAR AIR-SEA INTERACTIONS 

Heat lost from polar regions together with the heat source of the 
tropics form the basis of the global general circulation. The flux of heat 
| into polar regions from lower latitudes results from the general meridional 
temperature gradient from polar regions toward the equator due to unequal 
latitudinal heating of the earth's surface. Most of the heat advected into 
polar regions by the atmosphere is lost to space by long wave radiation, 
thus forming a basic mechanism for maintaining planetary radiative 
equilibrium. 

Heat is also advected into polar regions by ocean currents, but due 
to the relative inefficiency of air-sea heat transfer and the relative 
Slowness of ocean currents, this transfer is roughly one-tenth that 
attributed to the atmosphere. 

Sea ice is a dominant feature of both arctic and antarctic seas; however, 
due to major differences in tne physical configuration of sea and land 
masses of these regions, it plays decidedly different roles. The arctic, 
being essentially an ice-covered ocean, is extremely sensitive to climatic 
perturbations since the areal extent and depth of sea ice is critical to the 
regional energy balance. A surplus of about one-third of the annual energy 
influx would remove the arctic ice pack [Maykut and Untersteiner 1969, 1971]. 
On the other hand, the antarctic continental mass, covered by 2000 to 4000 
meters of ice, is relatively insensitive to small energy flux perturbations, 


although its sea-ice cover would be affected similarly to that of the arctic. 








Of major import to polar air-sea interactions is the degree to which 
the sea surface is covered by ice. Sea ice acts analogously to a ‘lid' 
on the sea, modifying the dynamics and continuity of the sea and the air- 
sea system. Sea-ice cover dramatically reduces the heat transfer between 
ocean and atmosphere by effectively limiting that transfer to a process of 
molecular conduction. Further, by intercepting the vertical flux of 
momentum transferred from the mean wind field to the sea surface by the 
Reynold's stress, it suppresses wind mixing. Finally, it reflects a 
considerable portion of the incoming short-wave radiation, imposes an 
upper limit on the sea's surface-boundary temperature, reduces latent heat 
transfer by impeding evaporation, and, through its changes cf phase, serves 
as a buffer, or damping mechanism. Alternatively, heat transfer processes 
are significantly enhanced over open-water areas in proximity to large ice 
masses, particularly during the winter, as surface air temperatures drop 
far below surface water temperatures. Wind and water stresses, in addition 
to directly affecting ocean-atmosphere heat transfer processes, keep the 
ice in nearly continuous motion, thus locally affecting heat transfer as 
well as producing mechanical Perce ene resulting in the formation of pressure 
ridges, polynya, and open leads. 

Polar air-sea interactions are therefore characterized by complex feed 
back wherein the ocean and atmosphere individually and collectively influence 
the thickness and areal extent of the ice, which in turn has a pronounced 


effect on the exchange in question. 


B. SURVEY OF PREVIOUS RESEARCH 
Range and depth of interest in arctic research has expanded considerably 


since Nansen's epic drift in the FRAM (1893-1896). From pioneering studies 








until only recently, research emphasis has been primarily of a descriptive 
nature, and it was not until the Soviet establishment of the drifting 
Station "North Pole 1 (NP-1)" in 1937 that a systematic study of the 
arctic interior was initiated. Launching of a comparable long-range 
arctic research program by the United States would await the establishment 
of "Fletcher's Ice Island (T-3)" in 1952. 

Attempts to understand the response of sea ice to environmental 
changes has long stimulated the interest of arctic researchers. Quantita- 
tive prediction of these responses has classically involved two primary 
approaches, the empirical and the theoretical. Although providing minimal 
insight to the dynamics of physical processes occurring between air, ice, 
and sea, empirical formulae relating sea-ice growth processes to observed 
Surface temperatures have accomplished basic predictive goals. The most 
sophisticated example of this approach is perhaps one advanced by Bilello 
[1961]. Theoretical studies involved analytical solutions and were there- 
fore severely limited in scope and complexity. One of the foremost of 
these analytical treatments is attributed to Kolesnikov [U.S. Naval Civil 
Engineering Laboratory 1966]. 

In the early 1950's the tenor of arctic research began to reflect 
increased emphasis on the theoretical approach to arctic dynamics as 
comprehensive data from systematic studies became available. It was at 
this point that studies relating the arctic energy balance to the mass 
budget and subsequent research directed toward heat transfer processes 
began to emerge. Development of physical relationships was confined to 
empirical and statistical approaches, or theoretical problems to which 


Simple, analytical, approximate solutions could be developed. 
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Classic arctic heat budget studies were presented by several authors. 
Badgley [1961, 1966] summarized and collated available data to deduce a 
typical heat budget at the surface of the central Arctic Ocean. Fletcher 
[1966] discussed annual patterns of atmospheric heat loss for each component 
of the heat budget in relation to the general atmospheric circulation, for 
both ice-covered oceans and ice-free oceans. Vowinckel [Arctic Meteorology 
Research Group 1964] independently investigated the amount of heat conducted 
through the Arctic Ocean ice to check previously calculated values for the 
heat released from the Arctic Ocean. Doronin [1966] presented an analysis 
of heat-balance components of the surface layer in the Arctic Ocean. 
Coachman [1966] proposed a model to explain the observed features of the 
oceanic regime of the Arctic surface layer during the winter. Muench 
[Baffin Bay - North Water Project 1971] discussed the physical oceanography 
of the northern Baffin Bay region in considerable detail. 

It was not until the 1960's and the attendant widened use of the high- 
Speed digital computer as a research tool that physical processes could be 
modeled free of constraints imposed by analytical solutions. Using finite 
difference schemes, a system of integral-differential equations could now 
be solved rapidly within acceptable error limitations. Typical of first- 
generation numerical models in the field is that developed by Untersteiner 
[1966], which predicted thickness and temperature of sea ice, thus allowing 
further insight into the relationship of climatic change to ice production, 
ice decay, and equilibrium ice thickness. The advanced one-dimensional 
thermodynamic model of central Arctic sea ice developed by Maykut and 
Untersteiner [1969, 1971] is a prime example of the present theoretical 


approach. 








Co PRESENT RESPARC Ea DURECTION 

A great deal of research effort is presently being directed toward 
developing numerical models of the atmosphere and of the ocean as well 
as atmosphere-ocean models, with the objective of ultimately producing 
an efficient, reliable global system model. Problems inherent in the 
realization of this objective are numerous and include not only a lack of 
refinement in respective models of the atmosphere and ocean, but also a 
lack of refinement in the technique of coupling the respective models. 
Further, observational input data is in general deficient for these models, 
particularly in the Arctic. Present state-of-the-art examples of models 
of this type are well represented by the sophisticated models of Manabe 
[1969] and Bryan [1969]. 

One of the more ambitious research projects in recent years is the 
Arctic Ice Dynamics Joint Experiment. AIDJEX is essentially a United 
States - Canadian cooperative effort to gain quantitative understanding 
of the interaction between atmospheric sea ice and fluid ocean fields of 
motion. fae ocKess at the time of this writing, AIDJEX is unique in its 
observational approach since, in order to gain a full understanding of 
interacting fields, it involves the collaboration of many research groups 
in making coordinated measurements of interacting fields of motion, stress, 
and strain - observing appropriate time and space scales - over a minimum 
period of one year. The fruits of this ambitious and imaginative study, | 
in terms of fresh insights, controlled data, and international cooperation 


appear to be vast already. 
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D. CONCEPTUALLY-RELATED STUDIES 

Studies of Arctic dynamics, energy balances and mass budgets range 
over a wide spectrum. Scientific literature reflects numerous experimental 
and theoretical models of microstructure and micro-scale energy transfer 
processes between ocean and atmosphere, ice and atmosphere, and ice and 
ocean, essentially occurring in boundary-layer regions. Pertinent to the 
work presented herein is research performed by Foster [1968, 1969] in which 
haline convection induced by the freezing of sea water was studied, and by 
Lake and Lewis [1970] who studied salt rejection by sea ice during growth. 
Work performed by Lewis and Walker [1970] illustrates seasonal cnanges in 
the temperature and salinity profiles beneath annual sea ice cover, and 
Stern and Turner [1969] reported on the forination of salt fingers and 
convecting layers. 

Directly related to the model proposed herein is the macroscale advanced 
model derived by Maykut and Untersteiner [1969, 1971] to forecast ice 
temperature and thickness, which defines the role of each component of the 
energy budget in relation to its effects upon the ice. Although comprehen- 
Sive in its treatment of physical processes related to components of the 
energy budget, the model admittedly has several shortcomings. Physical 
limitations are due to uncertainties in environmental data and include 
the necessary assumption of a constant oceanic heat flux and generally 
course albedo and energy-flux data. Limitations, arising from an inability 
to model certain processes include: 1) the neglect of mechanical stresses 
produced by winds and surface currents, 2) treatment of turbulent fluxes 
at boundaries as being independent of growth rates or the physical states 


of ice and snow, 3) imprecise treatment of downward heat transport from 
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melting snow, 4) undefined shapes of salinity profiles for various ice 
conditions, and 5) neglect of heat storage by melt ponds. Leads were 
considered only so far as they indirectly affected the energy fluxes over 
the bare ice. In spite of its limitations, this will certainly be a 
cornerstone for future thermodynamic ice models and provides a number of 


conceptual bridges used in the present work. 


BemeHESTIS OBJECTIVE 

It is Known that significant heat transfer occurs between the atmosphere 
and the ocean through open leads in the Arctic sea ice cover. As found 
experimentally by Badgley [1961], sensible neat toss is at least two orders 
of magnitude greater from open leads than from perennial sea ice. 

With the need for further refinement in models of the ocean and the 
atmosphere, as well as in model coupling schemes, the requirement for 
mesoscale examinaticn of ocean-ice-atmospnere neat-itransfer processes 75 
immediate and apparent. With appropriate observational input data becoming 
increasingly. more abundant, informative numerical models of mesoscale 
processes appear promising for furthering the understanding of physical 
processes and for realizing control of mass-energy fluxes in the arctic. 

The objective of this study was to develop a thermodynamic model of an 
open lead in central arctic sea-ice in which hydrodynamic processes of 
convection, advection, and diffusion of heat and salt were represented. 

By integrating the model numerically, the effects of various representative 
advective parametizations, and temperature and salinity profiles in the 
ocean's surface layer, on the heat flux through the lead and on the nature 
of ice formation, could be investigated as the lead reclosed by thermal 


processes alone. 
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TT. NATURE OF THE PROBLEM 


Since sea ice essentially insulates the ocean from the atmosphere and 
Suppresses important thermodynamic processes, it, along with heat-advecting 
ocean currents, becomes a primary factor in determining the vertical 
Structure in Arctic seas. However, the dynamic transition from ice-free to 
ice-covered seas is not fully understood. In particular, the nature of 
the fluxes of salt and heat, as well as the structure and depth of unstable 
convective overturn during initial phases of the freezing process, require 


Study. 


A. DYNAMICS OF THE SEA WATER FREEZING PROCESS 

The structure of the active surface layer of the Arctic Ocean is 
generally horizontally uniform with a very stable density gradient. How- 
ever, during the freezing process the density gradient may be significantly 
altered by three primary mechanism: 1) progressive cooling of surface 
waters from eeave 2) brine expulsion by the ice during the freezing pro- 
cess, and 3) advection of heat and/or salt by surface currents. 

The fundamental physical relation which makes the sea water freezing 
process unique may be understood by examining the curves of maximum density 
and freezing point for sea water, plotted with respect to temperature versus 
chlorinity or salinity (Figure 1). Specifically, for waters of salinity 
greater than 24.70 parts per thousand, the temperature of maximum density of 
the liquid state will always be lower than the liquid's freezing-point 
temperature under pressures existing in nature. Any process that tends to 


produce a vertical density instability will induce convective overturn in 


ifs: 





Figure |. 
Temperatures of Sea Water 
Maximum Density and 
Freezing-Point 

VS 
Chlorinity 


*after Sverdrup, Jonnson and Fleming [1942] 





Maximum Density 


0 
-| 
TEMP 
Freezing-Point 
O 
C 
-2 
-3 
-4 





3 10 iis 20 


Chlorinity 0/00 





the water column, resulting in mixing of the unstable segment such as 

would be produced in the freezing process by surface cooling and freezing- 
associated brine expulsion. Convective overturn will therefore be continuous 
in the unstable segment of the liquid water column from the onset of surface 
cooling to the point of total change-of-phase of the water column. A some- 
what less obvious but straightforward point should be mentioned in qualifi- 
cation of the foregoing. From an arbitrary water temperature above the 
freezing point, down to, but not including the freezing point, convection 

is temperature controlled during cooling since there is no ice formed, 
hence, no brine expulsion. During cooling of surface waters at the freezing 
point for existing surface salinities, convection is primarily salinity 
controlled since brine is being continuously expelled from newly forming 
ice. Within this latter water temperature regime, temperature changes play 
a minor role in effecting changes in surface water densities. 

As convective overturn is a mechanism for inhibiting freezing by 
continuously transferring heat from depth to the active boundary at which 
freezing is taking place, so it is also a mechanism for effectively removing 
heat from depths significantly greater than could be effected by molecular 
diffusion. An equilibrium ice thickness is reached when the resultant heat 
flux upward through convection, advection and diffusion in the water column 
balances the upward resultant heat flux through the ice, at the ice-water 
boundary. Prior to this point, ice is accreted in deptn since the upward 
heat flux through the ice exceeds the upward heat flux in the water column 


arriving at the ice-water boundary. 
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B. DISTRIBUTION OF VARIABLES 

The distribution of temperature, Salinity and velocity are considered 
within a volume of water under a refreezing lead. The water is assumed to 
be inviscid, horizontally homogeneous, horizontally isotropic, irrotational , 
and incompressible. Accelerations resulting from the coriolis effect and 
horizontal pressure gradients are assumed to be negligible. Momentum, heat 
energy, and salt mass are conserved within the volume such that no velocity, 
temperature, or salinity sources or sinks exist except at the volume's 
boundaries. Furthermore, currents are assumed to be steady-state and 
uniform both horizontally and vertically. 

Using the classic box-model approach, these conservative fields are 
represented continuously within the volume by individual continuity 
equations. For the vector field of velocity, where K is the diffusivity 


and § 


Sy 1S a source tern, 


= e + 
DW/ ap =H VKYVt S, Cl) 
which reduces to. 
v(t,x,y,z) = V = constant (2) 


i.e., the velocity is specified as uniform and steady. 
Similarly, the scalar fields of heat and salt, here collectively 


depicted by the scalar quantity @, are represented by 


DM/pe = Ve KVOTS, (3) 


where K is the diffusivity and 7 is a source term. 


3: 





This reduces to 


30, = -vevb+ Kw D ty . K(z) v9 (4) 


where Ku is the horizontal diffusivity and Vy is the vertical del operator 
Since, 

K(x) = K(y) = Ku = constant 

3g = 0 
The local, or time-rate-of-change, is therefore equal to the sum of the field- 


rate-of-change and the diffusive-rate-of-change. 


C. HYDRODYNAMIC PROCESSES 
Mechanisms for the hydrodynamic transport of heat and salt within the 
lead are advection and diffusion. However, from the dynamics of the sea- 
water freezing process, it is Known that convection plays a major role in 
transporting ooth neat and sait in the water column. Although it is not 
represented by a discrete term in the foregoing continuity equation, the 
convective transport term may be thought of as simply a temporary (when 
vertical density instabilities arise) modification of the vertical diffusion 
term taking the form 
my 8h. kal ae (5) 
K(z)+ 0 


which in essence results in free vertical (convective) transport. 


D. MASS AND ENERGY BALANCE 

With the foregoing analyses of continuity relations and mechanisms of 
hydrodynamic transport existing within the box meael of the open lead, the 
mass and energy balances existing at the volume's boundaries may be 


examined. 


Us, 





Over an arbitrary surface area of a reclosing lead, a given change of 
heat, 4Q, at the water surface must equal the heat added by advection and 
diffusion at upstream boundaries, 0 plus the heat added from depth in the 
lead by convection, advection, and diffusion, Q,, plus the heat added by 
surface ice accretion, Q¢, plus the open-water heat loss at the surface by 
conduction, evaporation and radiation, Q. anes plus the heat loss through 
the ice cover by conduction, Q. minus the heat diffused across vertical 
boundaries normal to advective transport, Q,> minus the heat transported 
across the downstream boundaries by advection and diffusion, Qqg, such that 

AQ = Qy + Qy + Qe + W(Q.tQ,40,) + Q;(1-W) - Q - Q (6) 
Where Wis the fraction of open water comprising the surface of the lead at 
any given time. The heat added from depth in the tead is heat ultimately 
derived from the Bering Sea Water (in the Canadian Basin) and from the 
Atlantic layer. 

By similar analysis, a given change of salinity, 4S, at the water's 
Surface must equal the salt added by advection ana diffusion at upstream 
boundaries, Si. plus the salt added from depth in the lead by convection, 
advection, and diffusion, Soe Plus sthe salt added by surface ice accretion 
(salt rejection), Se, plus the salt lost through surface ice melting (salt 


dilution), Se plus the salt added by surface evaporation, S,, minus the 


e 
salt diffused across vertical boundaries normal to advective transport, 


S,, minus the salt transported across downstream boundaries by advection and 


n 


Gigtusion, Sq» such that 


SS SS aaseey (1-W)(S_tS__) + W(S,) = Stages | (7) 


where again W is the fraction of open water comprising the surface of the 
lead at any given time. The salt added from depth in the lead is derived 
from the generally strong positive halocline existing at depths below 25 to 
50 meters. 
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Qualitatively then, before the equilibrium sea ice cover has been 
reestablished, there will exist differentials in the temperature and 
Salinity profiles of water columns. A primary process of interest in the 
refreezing lead is that of surface ice formation, which lends a unique 
character to lead hydrodynamics. The heat added/lost by surface ice 
accretion/melting is a function of the incremental change in the mass of 
ice, Al, and the heat of formation of sea ice, L., such that 

Qe = -Qn = LyAl | (8) 

Coupled directly to this process is the process of brine expulsion (salt 
rejection) by the sea ice during growth. The salt added/lost by surface ice 
accretion/melting is a function of the incremental change in the mass of ice, 
Pye che salinity of the sea ice, S;, (which in turn 1s dependent uponsrate- 
of-freezing, ice temperature, and age of the ice), and the salinity of the 
parent sea water, 5S,,, such that 

Sa SS SAS oo Sh) (9) 
Qualitative analysis of the salt rejection process by constituent salts is 
not appropos to either the scale or the nature of the study at hand; there- 
fore salt rejection was treated simply as a process of water-salinity 


modification. 
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III. EXPERIMENTAL PROCEDURE 


A. FORMULATION OF THE MODEL 
]. Establishment of a Representative Cross Section 

In deciding upon what physical dimensions would be appropriate for 
the model, it appeared most suitable to represent a small open lead in 
central arctic sea ice such as one would expect to observe between early 
to late winter. Therefore, the width of the lead was chosen to be 150 
meters, with an additional 50 meters under equilibrium sea-ice cover on 
the downstream side of the model lead to facilitate studying resultant 
thermodynamics and hydrodynamics in that region. 

The model's depth extent was set at 50 meters based upon lower 
boundary condition considerations. Over much of the Arctic Basin the 
vertical temperature gradient at that deptn is zero or very slightly posi- 
tive. Since a no-flux lower boundary condition was desirable in order to 
better study the effects of advection of conservative constituents in the 
lead, that depth was chosen. 

It was felt that the model could be simplified without the loss of 
any essential physics if derivatives along the lead (Y-axis) were set to 
zero. This reduced the problem to one of looking at a cross-sectional 
distribution in the lead. Figure 2 is a pictorial depiction of the open 
lead. 

2. Adaptation of the Continuity Equations to the Model 
Using the simplification of the previous paragraph which allows the 


use of a cross-sectional area rather than a true volume to represent the 
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lead, the horizontal advection and horizontal diffusion terms describing 
changes along the Y-axis were eliminated from equation (4). 

The problem was further simplified by assuming that all horizontal 
diffusion was negligible relative to horizontal advection, and that verti- 
cal advection could be disregarded. 

Applying the foregoing simplifications to the expanded form of 
equation (4), and choosing the state variables of temperature and salinity 
as the relevant dependent variables, the final form of the equations used 
to describe the conservative scalar fields of temperature and salinity 


within the lead was 


3A jy, = - UD/gx * Vaz K(E) BB/y, (10) 


where, 
P= 9 (t,x,2) 
3. Initial Conditions 
Sea-water temperatures and salinities within the lead at the 
instant that the lead is generated are represented by local stable tem- 
perature and salinity profiles previously existing beneati the equilibrium 
Sea-ice cover. The temperature and salinity fields are both horizontally 
uniform and vary only with depth in the lead such that 
T(0,X,Z) = T(Z) (11) 
Se. 2) = 252) , (12) 
The surface of the newly-generated lead is initially taken to be 
ice-free over the left (upstream) 150 meters, and ice-covered with equili- 
grium thickness ice over the right (down-stream) 50 meters. Surface 
temperatures are initially at the freezing point, in temperature equilibrium 


with the preexisting sea-ice cover. 
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4. Boundary Conditions 


The model's left, or upstream boundary, may be regarded as a 
source boundary across which sea water having a given temperature and 
salinity profile flows into the lead. Sea-water temperature and salinity 
profiles are taken as steady-state functions of depth alone, identical to 
the lead's initial temperature and salinity profiles such that 


T(t ,0,Z) 


tt 
| 
~— 
IN 
~~ 


(13) 
S052 


| 
—~/ 
-— 
IN 
~~ 


(14) 
The model's surface boundary is the primary region of mass and 
energy fluxes and is therefore the most dynamically active boundary. The 
boundary condition on temperature across the entire surface is that of an 
upward (or downward) diffusive flux of heat across the water surface such 


that 
O/- | T(t,X,0) = F(AT.2Z;,S; Ky Ky ) oS) 


where AT is the air-water temperature difference (an externally specified 
parameter), Plo eeMe moll idee Ge seme mess. S. is the surface ice salinity 
(an externally specified constant parameter), and K, and Kj are the co- 
efficients of thermal eddy diffusivity of sea water and thermal conducti- 
vity of sea ice respectively. When no ice is present at the surface, the 
Surface heat flux is taken as the cumulative open-water sensible, latent, 
and radiative heat fluxes, which are specified in the model as external 
parameters according to the air-water temperature difference. Since the 
model's surface at all times represents the water surface, the fluxes at 
that boundary are matched as ice is formed, such that the flux through the 


water surface 1S always equal to the flux through the ice to the atmosphere, 


but at no time does it exceed the maximum open-water flux to the atmosphere. 
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The boundary condition on sea water salinity at the surface 
specifies a no-diffusive flux condition between air and water as well as 


between ice and water and ice and air such that 


O/, S(t,x,0) = 0 (16) 


However, salt rejection from ice to water during the freezing process (as 
well as salt dilution during melting) is represented as a discrete process 
of physically adding the change in salinity produced by sea-ice accretion 


or melting to the surface waters such that 


Site. 0) = Sluts: t(AalI,4S) (17) 
where AI is the incremental change in the mass of sea ice, and AS is 
the salinity differential existing between the sea ice and its parent sea 
water. 
The bottom boundary, as previousiy discussed, is characterized 


as a no-flux boundary, thus 
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0 (19) 


Of, S(t.X,h) 
where h is the model's depth extent. 

The right, or downstream boundary, is also characterized as a 
no-diffusive flux boundary where whatever temperature-salinity field is 
generated by the thermodynamic processes of the lead is simply advected 
out. However, as mentioned previously, the right 50 meters of the lead 
model's surface is overlain by an equilibrium sea-ice cover. This boundary 
condition is therefore representative of a return towards the equilibrium 


vertical profiles. At the model's right boundary then 
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where W is the model's width. 

Figure 3 describes processes considered in the model cross section. 
Figure 4 schematically represents the model cross section with its governing 
continuity equations and boundary conditions. In this figure, the surface 
boundary condition is simply a statement of the balance of fluxes at the 
ice-water boundary. With no ice cover (6=0), the upward heat flux through 
the water surface, FS is equal to the sum of open-water sensible (pe). latent 


(F.), and radiative ue) heat fluxes. When an ice cover is present (670), the 


e 
upward heat flux through the water surface (2 =0t) is equal to the upward heat 
flux through the ice (2 =0 ) minus the heat added by surface ice accretion. 
The above relations are presented in terms of temperature gradients in Figure 
4 as they are applied to the governing continuity equation for temperature. 
Symbols not previously defined are p> volume density of sea ice; Gas heat of 
formation of sea ice; 6 ice thickness (volume of ice per unit surface area). 
5. The Finite-Difference Scheme 

A purely explicit finite difference scheme was selected to numeri- 
cally represent the linear parabolic partial differential equation (10). 
The time rate-of-change term was represented by forward differencing, the 
field rate-of-change term by backward (upstream) differencing, and the 
diffusive rate-of change term by central differencing. Thus over the interior 


(non-boundary) region of the model 


n+] n n n 
_ gf n I -29. +0. 

Dik Ds ik — Dek - Di] ,K +k Paka Pe ih Ioko (22) 
At 4 Xx : (4z)2 


?7 





Figure 3. 


Description of Processes Considered in Mode] 





(1) 


(2) 


(3) 


(4) 
(5) 


Advection of steady-state temperature-salinity field into the 
lead from under equilivrium ice pack. 


Heat loss to the atmosphere as a function of ice thickness, 
ice accretion/melting, and salt rejection/dilution. 


Heat and salt transport by convective overturn, advection and 
eddy diffusion. 


No transport across boundary. 


Advection of temperature-salinity field generated by lead 
thermodynamics out of the lead model. 
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where, 


n a, eo 
Pek JAX, KAZ (23) 


At the left boundary, equation (22) was inapplicable since tempera- 
ture and salinity were both steady-state functions of depth alone. There- 


fore, at this boundary at all times 


(24) 


At the surface boundary, the diffusive term for temperature was 


modified in equation (22) since 


nN nN 
V5 ket 7 5k 


= VT (25) 
242 
or, 
i =e era (a (26) 
ok-1 O° 5,kH a a 
The diffusive salt term at the surface was modified in equation (22) 
since ; 
n n 
Ce - Sa 
Hog] ip Gelleeuer (27) 
PLETE 
oar 
n oan 
Sjk-1 = $5 ,k4 (28) 


At the bottom boundary, the diffusive rate-of-change terms of both 


temperature and salinity were modified in equation (22) since 


n n 
ah - Q. 
Pi sk# TG sk-1 - (29) 
ZAZ 
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Or, 
n Pahl 
Bs k+1 = 95 jk-1 ei) 


Finally, at the right boundary, equation (22) remains unmodified 
Since advection is the only transport mechanism operating across it. 
Stability analysis of the finite difference equation using the 
Fourier series method [Smith 1965] shows that 
n+] _ on U K 
Sue” ak [2 Hae * 2 wep t (31) 
where e€ represents an arbitrary error introduced into the problem. There- 


fore, the finite-difference equation (22) is stable when the absolute value 


of the amplification factor remains less than or equal to unity 





ac{o+2 | <1 (32) 
| LAX ez)e) 
See nepresenitat ion of Physical Processes | 

Convective overturn is treated nominally through an artifice of the 
vertical diffusion term as in equation (5). Numerically however, convective 
overturn and its associated mixing iS accomplished through arithmetic 
averaging of temperatures and salinities within regions of vertical density 
instability as determined by sigma-t values, followed by recomputation of 
the resultant sigma-t value for the newly mixed (and neutrally stable) 
region. Subroutine OVRTRN accomplishes convective overturn in the model. 

Advection of heat and salt (temperature and salinity) in the lead is 
driven by a current, u, in the positive-x direction. Values and depth pro- 


files of the current as used in the model are given in subparagraph B of this 


section. 
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Diffusion of heat and salt in the vertical are determined by the 
coefficients of vertical temperature and salinity diffusivities, both set 
at constant values of 10 cm* sec ~! in the model. 

Formation of sea ice is accomplished as in equation (8). More 
specifically however, when a given mass of water,4w, is cooled AT degrees 


centigrade below its salinity dependent freezing-point temperature, the 


incremental mass of ice formed, AI, 1s 
Al= “p(AWAT)/ (33) 


where Cy is the specific heat of sea water and L; is the heat of formation 
of sea ice. 

Salt rejection (or dilution) resulting from surface sea ice accretion 
(or melting) is described in equation (17). Specificaliy, if AI is the 
incrementa! change in the mass of sea ice and AS is the salinity differential 
existing between a unit mass of parent sea water and a unit mass of its 
Generated sea ice, then the incremental salinity change, SSp, in the surface 


waters resulting from that freezing or melting process is 


7. Physical Quantities and Relations 
Sea water chlorinity is related to salinity [Knudsen 1902] by the 


relation 
Cl=(s - 0.030)/1.8050 (25) 
Freezing-point depression is empirically related to sea water 
chlorinity by Thompson [Sverdrup, Johnson and Fleming 1942] by the relation 


Te, = -0.0966°C1 - 0.0000052:C13 (36) 





The density of sea water as a function of salinity and temperature 


(sigma-t) [Knudsen 1902] is empirically computed through 


¢ = #,(C1) (37) 
A, = f(T) (38) 
Br = f(T) (39) 
zy = €4(T) (40) 
a. = fe(Z7, 0), Az, Br) (41) 


The specific heat at constant pressure of sea water 1S empirically 
related to its salinity at normal temperature and pressure by Kuwahara 


[Sverdrup, Johnson and Fleming 1942] by the relation 


C, = 1.005 - 0.004136-S + 0.0001098-S¢ 
- 0.000001324-s3 (42) 
The thermal conductivity of sea ice was determined using the relation- 
Ship established by Untersteiner [Maykut and Untersteiner 1969, 1971]. 


7 _ BS(z) 
kK. - Ky 4 7 » 1<G (43) 
where, . 


the subscript f refers to pure ice, 

ke = 0.00486 cal/cm-sec-°C, 

B = 0.28 cal-cm“/g-sec, 

T = temperature of ice in degrees Celsius, 

S(z) = ice salinity at depth z in g/cm. 

The heat of formation of sea ice was represented by the simplified 


relationship [Pounder 1965] 


by = heals Soe) ma 
where, 

ey oe edly Ce 

o = Salinity of sea ice, g/kg, 


S = salinity of parent sea water, g/kg. 
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B. ENVIRONMENTAL DATA INPUTS 

Environmental data used in the model fell into two basic categories: 
that which remained constant throughout all runs of the model, and that 
which was variable from run-to-run. Representative of the former were 
initial equilibrium-ice thickness, set at 200 cm, an ice salinity of 8.0 
g/kg, and an ice density of 0.91 g/cm>, Further, representative open water 
sensible, latent, and radiative heat flux densities as presented by Muencn 
[Baffin Bay-North Water Project 1971] were applied in the model. However, 
as ice formed at the surface, the latent heat flux density was modified so 
as to decrease exponentially with increasing ice thickness. The radiative 
heat flux density was modified so as to decrease from open water values as 
ice formed, approaching a constant value representative of the radiative 
flux density over pack ice as also presented by Muench. The sensible heat 
flux density was modified internally by the relation shown in Figure 4 as 
surface ice thicknesses increased. All surface heat flux densities were 
representative of an air-water temperature difference of 25°C with no wind. 
Finally, the coefficients of vertical temperature and salinity eddy diffusi- 
vity were set at 10.0 cmé/ Sec. 

Primary variable parameters in the model were horizontal current magni- 
tudes, vertical temperature profiles, and vertical salinity profiles. 
Currents applied in the model were within the range of zero to 10.0 cm/sec. 
Extreme test case values were 0.0 and 10.0 cm/sec. Additional test case 
values of 2.0 and 4.0 cm/sec were used to study one particular temperature- 
salinity field combination with its associated hydrodynamics in further 
detail. A value of 7.0 cm/sec was chosen to represent a typical horizontal 
current occurring in the lead, and was subsequently used as the real case 


value. 
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Test case values of vertical temperature and salinity profiles 
representative of extremes in vertical gradient were applied. Profiles 
A represented strong gradients, while profiles B represented weak 
gradients (Table I). Real case temperature and salinity profiles, 
qualitatively representative of early winter (profile C) and late 
winter (profile D) conditions in the central arctic [Coachman and 
Barnes 1961, 1962], were chosen to represent the typical lead. 

Table I summarizes the profile designations of the three variable 


parameters and their various values. 


C. MODEL INTEGRATION PLAN 

In order to gain maximum insight to physical processes occurring in 
the lead, as well as to study the effects of current magnitude, vertical 
temperature gradient, and vertical salinity gradient on lead thermodynamics, 


the series of test variables was run in the model. Represented by run 


A) 


one through eight, these test cases involved running all possible combina- 
tions of current magnitude (u) profiles A and B, temperature (T) profiles 
A and B, and salinity (S) profiles A and B. These test cases were 
integrated over a lead reclosure period of two days. 

As a result of the interesting dynamics occurring in run 6, additional 
runs 9, 10, and 14 were scheduled using temperature and salinity profiles 
identical to run 6 (and 2), but with current magnitude profiles C, D, and 
F respectively, in order to further investigate these phenomena. Integra- 
tion was again for a period of two days. 

Finally, the real-case variables, representative of naturally 
occurring parameters, were put into the model and integrated aver two-day 


periods in runs 11 and 12, and over a 30 day period in run 13. 
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Table I. 


Environmental Variable Profiles 





Water Temperatures 


od VT, 


jH—_A__| 6 __} _¢_}_p __ 
0 -1.50 1.50 -1.68 ee2 
5 -1.40 -1.49 -1.68 = 1532 
10 =1 30 = 148 -1.68 -1.82 
15 5/640 =| ay -1.67 -1.82 
20 -1.10 -1.46 = GG =|. 
25 -1.00 1.45 -1.65 S| 86 
30 -0.90 eie4aa -1.64 -1.79 
35 -0.80 =e -1.63 =| 78 
me 49 | “—C/00)) See) | Sto | ae 


45 a eas: 
90 ; -1.795 





Water Salinities (g/kg) 


Depth Profile 
@) . lS aS aa 
0 Zoa00 28.00 1.00 5380! 
5 Cano 28.01 31.01 S612 
10 29.00 tae Ma S02 Sie 
Is 29.90 Bele: 3.206 35855 
20 30.00 28.04 lis Sou) 
25 Sl a0, Deeley heal 33.60 
30 31.00 28.06 See 35 85.03 
55 S]| 20 28.07 31.45 33.66 
40 S200 28.08 Soe) Boog 
45 62.30 20.09 S105 Bon?) 2 
50 Sse U0 io Mai Ble fo BOD 
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Table II presents a summary of runs conducted with their respective 


variable combinations. 


D. MODEL LIMITATIONS 

In addition to the assumptions made earlier in arriving at the 
simplified model of an open lead, basic limitations in the model's 
capabilities are: 1) the atmosphere is treated as an external parameter 
such that it affects lead dynamics but cannot be affected in turn by the 
lead; 2) sea ice is treated as a quasiexternal parameter since it is 
depicted as a simple, uniform mass floating on the water's surface. 
Density and salinity are constant throughout this mass in both space and time, 
and the vertical temperature gradient is taken to be linear throughout. Its 
thermal conductivity is determined only at the ice-water boundary , which is 


warmest and therefore limiting through the relation given in equation (43), 


“p 
Pir 


he 


(hb 
=4) 


applied to the upward heat-fiux density tnrougn 


" 


and an upper limit 4 
ace such that at no time is the flux through the ice allowed to exceed that 
over open water; 3) mechanical lead closure (such as by wind-drifted ice) 
is disallowed; 4) snow cover over the lead's surface is not modelled; 5) 
super-cooling of surface waters [Coachman 1966] is not portrayed; 6) as 
convection-induced mixing occurs instantaneously in the model, so also 

does surface ice melting whenever heat is brought to the surface by 
convection. This results in an artificiality in that, although heat and 
salt are conserved in the model, heat brought to the surface under an ice 
cover immediately melts ice until no further surface heat is available for 
melting. Therefore, under this limitation, it is not possible for surface 
heat to be advected downstream under an ice cover to promote downstream 


melting. Only in cases when no ice is present to be melted may this heat 
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Table I1. Schedule of Runs and Variaole Combinations 


Variable Profile Period of Integration 
Run U (days) 


Se ees 
1 A A A 2 
2 A A B 2 
3 A B A 2 
4 A B B 2 
5 B A A 2 
6 B A B 2 
7 B B A 2 
8 B B B 2 
9 i A B 2 
10 D A B 2 
i E ¢ C 2 
12 E D D 2 
13 e | 2 D | 30 
14 rp 6[ OA B 2 





be advected downstream; 7) in situations wherein vertical salinity gradients 
existing in the water column are weaker than a certain minimum gradient, 
the model will over-mix the water column in depth. Specifically, the model 
forms ice in response to surface cooling, thereby rejecting Salt into the 
surface water layer of depth Z. The model then mixes the induced un- 
stable region as determined by relative sigma-t values of density. If 
however, there is heat available at some shallower depth than that to 
which the water column was mixed which could have decreased the initial 
amount of ice formed, then although salt, ice, and heat are ultimately 
conserved, potential energy is not conserved. Initial freezing and resul- 
tant salt rejection have induced convection to a greater depth than would 
have actually been the case in a continuous process. Therefore, for 
vertical salinity gradients weaker than a certain minimum gradient, the 
model's numerical convective process acts as a potential energy source. 
Quantitatively, the limiting minimum vertical salinity gradient which can 


be accurately treated by this model, in g/kg per vertical grid space, 1s 


; > (See) = A1/ Az (45) 


where p. is the ice density, S,. is the salinity of sea water, 5. is the 


W 
salinity of sea ice, AI is the incremental change in the volume of ice, 
and AZ is the vertical grid spacing. In this model, the limiting gradient 
is in the neighborhood of 1 X 107° g/kg per vertical grid space, which is 


lower by an order of magnitude than the weakest salinity gradient used. 


So 





IV. PRESENTATION OF DATA 


Some results of integrating the model using the combinations of 
environmental variables summarized in sections III.B. and III.C. are 
presented in tabular and graphic form. Three primary output variables are 
treated: surface ice thickness, surface heat loss to the atmosphere, and 
depth of penetration of surface-freezing induced convective overturn. 
Throughout the section, these three output variabies are related by run 
number (section III.C.) to the input variables of current magnitude, 


vertical temperature profile and vertical salinity profile (section II1.B.). 


femeeS! CASES 

Ice thicknesses are treated in two ways: by time growth rate and by 
Spatial profile across the jiead. Values for growth rates were arrivec at 
by tabulating the maximum ice thickness occurring anywhere in the lead as 
a function of time. Values for spatial profiles were taken as the ice 
thicknesses across the model's surface (200 meters) after 48 hours of 
freezing. Values over the right 50 meters of the model's surface, initially 
covered by 200 cm of ice, are tabulated to show the downstream surface 
effects of lead thermodynamics on the equilibrium ice cover. 

Surface heat losses to the atmosphere are presented as cumulative heat- 
loss densities (heat loss per unit surface area) occurring at each grid 
point over the surface of the initially open lead after 48 hours of lead 


reclosure. Additionally, these surface heat-loss densities are averaged and 


tabulated separately as mean surface heat-loss densities. 
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Finally, maximum depths of penetration of surface-freezing induced 
convective overturn are tabulated. Due to the numerical nature of the 
model however, these depths should be treated as bounding depths rather 


than precise depths as would be indicated for continuous schemes. 


B. REAL CASES 

Ice thicknesses, surface heat losses to the atmosphere, and depths of 
penetration of surface-freezing induced convective overturn are treated in 
an analogous manner. Additionally however, ice tnickness and mean surface 
heat-loss densities are presented as functions of time for run number Is, 
in which the model was integrated to represent 30 days of lead reclosure. 
Values for surface ice thicknesses were determined analogously to those of 


48-hour ice-growth rates. 
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Surface Heat-Loss Densities, Test Cases 


Table V. 


48-Hour Surface Heat Losses (x 10° ACHE) 
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Table VI. Mean Surface Heat-Loss Densities, Test Cases. 


48-Hour Mean-Surface Heat Losses 
(x 102 cal/cm¢) 


Run 
Number 6 









OT! BPO WO WO MS © WO 


Om AHDOLP WN — 





Table VII. Waximum Depth of Penetration of Convective Overturn, 
Test Cases. 


Run Depth 
Number (m) 
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fabie 1X. Ice Thickness Profiles, Real Cases. 


Ice Thicknesses (cm) Across Lead Surface 
48 Hours After Lead Generation. 


Run Number 
m 11 2 


10 24.7 24.6 
20 24.7 24.6 
30 24.7 24.6 
40 24.7 24.6 
50 24.6 24.3 
60 24.3 24.0 
70 24.4 ZA 2 
80 24.4 24.2 
90 24.4 24.1] 
100 24.4 Zoo 
110 24.3 Figo) 
120 24.2 2Gt Gan \ i 
0 
0 
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Table X. Surface Heat-Loss Densities, Real Cases. 


48-Hour surtage Heat Losses 
(x 10¢ cal/cm°) 























X Run Number 
m) 12 
10 13.6 138 
20 L336 13.9 
30 136 1322 
40 [3.16 13.9 
50 Wor. 14.3 
60 14.2 14.7 
70 14.0 14.4 
80 ime 14.3 
90 ie6 455 
100 ees: 14.3 
1tG (icires: Pes 

133s 14.3 

eee 14.3 

lies 14.2 

1.99 Le 

1.99 

1.99 

1.99 

1.99 

1.99 
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Table XI. Mean Surface Heat-Loss Densities, Real Cases. 


48-Hour Mean Surface Heat Losses (x 10¢ cal/cm¢) 
Run ss 
Number AQ 


1] Tad 
12 14.2 


Table XII. Maximum Depth of Penetration of Convective Overturn, 
Real Cases. 


Run 
Number 














Depth 
(m. ) 





i 
ifs: 
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Table XIII. Ice Growth Rate, Run 13. 


Maximum Ice Thickness (cm) Occurring in the 
Lead's Surface at Five-Day Intervals. 


Ice 
Day ' Thickness 





Table XIV. Mean Surface Heat-Loss Densities, Run ee 


- 


Mean Surface Heat-Loss Densities (x10¢ cal/cm*) 
at Five-Day Intervals. 
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¥.  CONCLUSTONS 


Through the preceding data and the time-dependent temperature-salinity- 
density fields generated in the lead model, numerous insights to the process 
of lead reclosure were afforded, as well as to the hydrodynamic validity of 
the model under certain conditions. 

Perhaps most pronounced was the very strong role played by the near- 
Surface vertical salinity gradient in maintaining water-column stability. 

A strong halocline near the surface tended to insulate deeper waters in the 
lead from the thermodynamic activity occurring at the surface. Surface 
waters whose densities had been increased by salt rejected by newly forming 
ice, was not able to mix vertically until density instabilities occurred in 
the surface region, thus inducing convective overturn. Since the presence 

of a strong surface halocline would require higher surface-salinity accumula- 
tions and hence permit more ice formation and associated salt rejection 
before this surface layer would become unstable, a strong halocline would 

be expected to inhibit convective overturn, botn in frequency of occurrence 
and depth of penetration in the lead. Added to this insulating effect is 

the fact that, as surface salinities increase, surface freezing temperatures 
correspondingly decrease such that freezing rates are proportionately lowered. 

The role played by the vertical temperature gradient in lead thermodynamics 
was strongly dependent upon the existing vertical salinity gradient. Since 
density changes in sea water at temperatures near the freezing point are 
almost totally a function of salinity changes, the insulating qualities of 


the thermocline are undefinable when treated independently of the halocline. 





However, in situations wherein weak vertical salinity gradients (poor in- 
Sulation) and strong vertical temperature gradients coexisted, convective 
overturn resulted in relatively large quantities of heat being transported 
to the surface. Freezing and associated salt rejection were thus inhibited 
and surface melting was induced. In this manner, the temperature stratifi- 
cation was superimposed on the primary control exerted by the vertical 
Salinity gradient. 

The effect of horizontal current magnitude on lead Stability appearea, 
as in the case of the thermocline, to be pean ary to the effect of the 
halocline. Strong currents will advect surface salt accumulations out of 
the lead relatively quickly, and will in general allow less time for lead 
thermodynamics to affect a given parcel of water advectively transiting 
the lead's width. This allows less change in the parcel's physical charac- 
teristics. Through these effects, the current also acts as an insulating 
agent. 

Ice growth rates were clearly dependent upon combinations of the above 
three variable parameters. Neglecting spacial ice profiles across the lead 
temporarily and examining only the maximum ice thicknesses occurring on the 
lead's surface as a function of time, the hichest growth rates were exhibited 
for cases where a strong current, a weak thermocline, and a weak halocline 
occurred together. Conversely, iowest growth rates resulted from the co- 
existence of a strong current, a strong thermocline, and a strong halocline. 
Very little difference however, was noted in the high growth-rate case 
when currents were weak instead of strong, nor in the low growth-rate case 
when currents were weak or when currents were weak and the halocline was 
weak. Analysis of these results showed that high growth rates were in all 
cases associated with weak vertical temperature gradients and low growth 


rates with strong vertical temperature gradients. Further, that current 
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strength was relatively unimportant in changing growth rates cane 
strong haloclines were present, or when weak temperature gradients together 
with weak salinity gradients were present. 

The effect of the near-surface vertical temperature gradient appears 
then to dominate in determining surface freezing rates. Strong gradients 
resulted in lower rates of ice formation than in the case of weak gradients. 
Vertical transport of heat from depth in the lead by both eddy diffusion and 
convective overturn is then modified by both current strength and the near- 
surface halocline acting as insulating agents. 

Ice thickness profiles across the lead's width varied considerably with 
different combinations of the three basic variable parameters. Uniform or 
nearly-uniform ice thicknesses across the lead were in all cases associated 
with strong near-surface haloclines. Variations occurred where weak halo- 
clines existed, growing larger when weak haloclines and strong thermocl ines 
occurred together. The greatest surface irregularity, both by ampiitude of 
variation as well as by spacial frequency of variation, was found when weak 
currents, weak haloclines, and strong thermoclines occurred together. 

Analysis of these cases indicated that the halocline was the dominant 
variable in determining uniformity of ice thickness across the lead through 
insulation of the surface from the heat at depth in the lead. When convective 
overturn occurred as a result of increases in surface salinity, surface freezing 
rates were altered at the point of convective overturn and produced surface 
irregularities. Since rejected salt is being continuously advected down- 
stream, there will exist horizontal salinity gradients which, when exceeding 
the magnitudes of vertical salinity gradients, will induce convective overturn 


and resultant surface melting at that point. Downstream of this HOMER 
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surface waters are diluted due to the upstream melting, thus experiencing 
increased freezing rates since their freezing point has effectively been 
raised. 

The near surface thermocline appeared to affect surface-ice uniformity 
through determining the amount of heat transported to the surface during 
convective overturn, and thus the amount of melting and dilution occurring. 
Finally, current strength appeared to act analogously to the halocline as an 
insulating agent, though less effectively. By flushing surface-salt accumu- 
lations out of the lead more rapidly, a strong current appears to inhibit 
convective overturn and hence variations in ice uniformity. By decreasing 
the amount of change a given parcel of water experiences over a given hori- 
zontal distance in the lead, a strong current appears to produce less 
intensive but more extensive lee effects from surface dynamics. 

Surface heat-loss densities varied across the lead's surface in response 
to ice thickness profiles. Since open-water sensible, latent, and radiative 
heat-flux densities used in this work totalled 2.00 x 107¢ cal/cm¢-sec as 
compared to a corresponding value over the 200 cm of ice of 1.15 x 1073 
cal/cm@-sec, resultant heat-loss densities were considerably larger over 
regions of the lead surface covered with relatively thinner ice. Uniformity 
of surface heat-loss densities over the lead's surface may then be directly 
related to the uniformity of the ice profile over the surface, while relative 
magnitudes of heat-loss densities may be inversely related to relative ice 
thicknesses. 

The extreme case of lead instability and non-uniformity (weak current, 
strong thermocline, and weak halocline) requires special comment. A high 
mean surface heat-loss density characterized this case, which was expected as 


a result of its highly non-uniform surface-ice profile. ‘“owever, the lack 
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of evidence of these highly non-uniform ice covers occurring in nature in the 
absence of mechanical deformation forces prompted further analysis. By 
studying model output fields of runs 2, 6, 9, 10, and 14 (particularly the 
horizontal density gradients generated in the lead), an upstream pressure 
force was found to exist in the vicinity of the source boundary which was 
capable of performing on the order of ten ergs of work against a unit mass 

of water entering the lead at the source boundary. Since a thermodynamical ly 
generated force of this type would completely alter the hydrodynamics 
occurring in the lead for low current magnitudes, the model appears to be non- 
representative of actual lead thermodynamics when using very low current 
magnitudes with strong thermoclines and weak haloclines. Stronger currents, 
on the other hand, had generally weaker horizontal density gradients associated 
with them and were also proportionately less affected by these pressure forces 
Since their capacity to do work against the horizontal pressure force in- 
creases as the square of their magnitude. 

Results of run 14 proved to be of particular interest through the in- 
Sights they afforded to the preeess of Convective overturn. Since no GuiryenE 
was present, there was no horizontal coupling between model grid points. As 
a result, although freezing progressed uniformly over tne surface of the lead, 
strong differential horizontal density gradients developed along the ‘source’ 
boundary as cooling progressed. Of prime interest however, was the graphic 
illustration of the role played by the current as a convective-overturn 
scaling factor. With zero current, a strong temperature gradient arid a 
weak salinity gradient, freezing-induced convective overturn penetrated to 
a depth of 35 meters roughly 18 hours after lead generation. Without the 


flusning action of a horizontal current or the insulating effect of a strong 
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halocline, the effects of surface thermodynamics were felt much deeper in 
the lead. It appears probable then, that jn cases wherein weak currents and 
weak haloclines exist in open leads (marginal vertical stability), surface 
freezing will induce, through convective overturn and vertical mixing, sub- 
stantial vertical and horizontal gradients in both temperature and salinity. 
As a result, significant fluxes of heat, salt, and momentum would be genera- 
ted in the waters of the lead, and sophistication beyond that possessed by 
the work at hand would be required to accurately represent the resultant 
lead hydrodynamics. 

Summarizing the analysis of results from integrating the open-lead mode | 
over the experimental ranges of variables, the near-surface vertical salinity 
gradient appears to exert the dominant force in maintaining vertical density 
stability in the water column through its insulating effect. By inhibiting 
convective overturn, a strong haloc!ine was associated with horizontal ly- 
uniform ice thickness profiles and consequently uniform surface heat-t lux 
densities. Current strength and the vertical temperature gradient appeared 
to be secondary as insulating agents. 

The near-surface temperature gradient appeared to be the dominant 
variable in controlling the rate of ice formation Since its strength 
directly influences the amount of heat transported ae jce-water boun- 
dary from depth by both eddy diffusion and convective overturn. In this 
case, the halocline and the current strength appeared secondary in their 
effects on the freezing rate. 

Mean surface heat-loss densities were higher for more non-uni form 
surface ice thicknesses, which in turn may be inversely related to the 


strength of the near-surface halocline. This conclusion could prove to be 
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of significance to critical examinations of the heat buaget over large 
areas of the arctic in which significant regional variations in the near- 
surface halocline might reasonably be expected. 

Representative freezing rates, ice-thickness profiles, and surface 
heat-loss densities found in integrating the model using real-case variables 
agreed well with values in the literature. Due to the similarity of current 
strengths, near-surface thermoclines, and near-surface haloclines, the fore- 
going relations were generally not discernable, although the effect of down- 
stream salt accumulation was noticed in slightly decreasing downstream 1ce 
thicknesses. 

Finally, it is concluded that, with the exception of cases in which 
weak horizontal currents occur together with strong thermoclines and weak 
haloclines, the model presented herein provides an accurate physical 
representation of the thermodynamics occurring in an open jead in central 


arctic sea ice over the experimental range of variables. 
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VI. RECOMMENDATIONS FOR FURTHER RESEARCH 


Although the surface boundary of the model could be considerably more 
sophisticated in its treatment of energy fluxes and transport processes, 
models already exist which have this sophistication, notably that of Maykut 
and Untersteiner [1969, 1971]. Coupling the model at hand to their model 
would therefore seem a more reasonable approach to refinement of the surface 
boundary in order to model air-ice-sea interactions and mass-energy balances. 

The hydrodynamics occurring in a refreezing lead however, are still not 
accurately modelled, nor are they fully understood despite the insights 
afforded by the research performed. Extended development of this work 
would promote further conceptual understanding of these processes as well 
as provide a more physical ly-accurate mode?. Incorporation of accelerations 
resulting from thermodynamically-generated horizontal pressure gradients, 
as well as the coriolis acceleration, would considerably refine the model. 
Additional refinement would be gained by removing previously discussed 
limitations on: 1) the minimum vertical salinity gradient able to be 
accurately treated, and 2) the preclusion of downstream advection of warmer- 
then-freezing-point waters under an ice cover. 

In summary then, further research could profitably be applied to the 
work at hand to further refine the hydrodynamic accuracy of the model, as 
well as to couple it to larger-scale models of air-ice-sea interactions and 


mass-energy balance predictive schemes. 
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APPENDIX A 
COMPUTER PROGRAI DESCRIPTION 


The computer program is written in FORTRAN IV and was used with the 
IBM 360/67 computer system at the W.R. Church Computer Center, Naval 
Postgraduate School. The program has been generalized such that sea 
water temperature and salinity profiles used as initial and left boundary 
fields are treated as input data. A design feature of the program allows 
the use of an internally variable time step to facilitate adjusting the 
model's temporal resolution. Additionally, the model allows variable grid 
Spacing as well as variable lead dimensions, requiring changes only in the 
DIMENSION, COMMON and FORMAT statements. Although all currents used in 
the model were constant in time and space, the model allows use of both a 
Geptn-dependent and time-dependent horizontal current, again internaily 
variable. With only minor modifications, the program could be easily ex- 
tended to include a bottom-boundary heat and salt flux, more sophisticated 
treatment of air-sea interactions at the surface boundary, and feedback of 
accelerations, resulting from horizonta! density gradients generated by lead 
thermodynamics, to the motion of water through the lead. 

The program is divided into two basic sections: the main program and 
Subroutine OVRTRN. The main program is further subdivided into seven parts: 
problem initialization, data input, computation of time-dependent temperature 
and salinity fields, ice-accretion/melting and salt rejection/dilution, with 
vertical stability check and computation of resultant heat flux densities, 
filling temperature and salinity arrays with resultant values, computation 


of applied surface heat flux densities and conversion to an integrated surface 
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temperature gradient, and writing output. Subroutine OVRIRK accomplishes 
vertical convective overturn and mixing of unstable regions in the water 
column. 
FORTRAN IV symbols for primary program parameters are defined below: 
AH Downstream-advected heat-flux density (cal/cm-sec) 
AHF Applied surface heat-flux density (cal/cm¢-sec) 
BTC Sea ice thermal conductivity constant (cal-cm*/g-sec) 
CAH Cumulative downstream-advected heat density (cal/om®) 
CAHF Cumulative applied surface heat density (cal/cmé) 
CHC Cumulative subsurface cooling heat density (cal/cm) 
CHF Cumulative applied surface heat density less cumulative 
| downstream-advected heat density (cal/cm¢) 
CHIF Cumulative heat density of surface-ice formation (cal/cm) 
Gi Thermal conductivity of sea ice (cal/cm-sec-°C) 


OL Sea-water chlorinity (g/kg) 


CO Thermal conductivity of pure ice (cal/cm-sec-9C) 
CP Specific heat of sea water (cal/g-°C) 
DI Rate of surface ice formation (cm3/cme-sec) 


DSALT Salt rejected per cm? of ice formed (g/kg) 
DTAU Time step (sec) 

DTEMP Air-water temperature difference (°C) 

DX Horizontal grid spacing (cm) 

DZ Vertical grid spacing (cm) 


2_sec) 


ees Allowable rate-of-freezing error (cm?/cm 
ERR Allowable ice thickness error (cm?/cm2) 


CRADI Vertical ice temperature gradient (°C/cm) 


7A 





H 
HC 
HF 


HIF 
HO 
HS 
MH] 
MV I 
MVHI 


QI 


QL 


QR 


QS 


RHO 
RHOI 
S 
SALI 
SI 
SIE 


Model depth extent (cm) 

Subsurface cooling heat-flux density (cal/cm¢-sec) 
Applied surface heat-flux density less downstream-advected 
heat flux density (cal/cm¢-sec) 

Heat-flux density of surface-ice formation (cal/cm¢-sec) 
Heat of formation of pure ice (cal/g) 

Heat of formation of sea ice (cal/g) 

Number of grid columns 

Number of grid rows 

Total number of grid points 

Heat-flux density through sea ice to the atmosphere 
(cal/cm¢-sec) 

Open-water latent heat-flux density to the atmosphere 
(cal/eme-sec) 

Open-water radiative heat-flux density to the atmosphere 


(cal /cm¢-sec) 


-Open-water sensible heat-flux density to the atmosphere 


(cal/cm@-sec) 

Density of sea water (g/cm”) 
Density of sea ice (g/cm) 
Sea-water salinity (g/kg) 
Sea-ice salinity (g/kg) 
Sea-ice salinity (g/cm?) 


Sigma-t 


SSCHG Salinity change in surface layer of depth DZ (g/kg) 


Ti 





SSFLUX Salinity change in one em? of sea water at the ice-water 
boundary (g/kg) 
STGD Total applied surface temperature gradient at the water's 


surface (9C/cm) 


if Sea-water temperature (°C) 
Tau Time (sec) 
TF Sea-water freezing temperature (°C) 


TGRADL Surface temperature gradient due to latent heat-flux 
density (9C/cm) 

TGRADR Surface temperature gradient due to radiative heat- flux 
density (9C/cm) 

TGRADS Surface temperature gradient due to sensible heat- flux 
density (9C/cm) 

it Temperature of lower ice surface (9C} 

U Horizontal current magnitude (cm/sec) 

VEC ‘Vertical eddy conductivity (cal/cm-sec-°C) 

VSK Vertical salt eddy diffusivity (cm/sec) 


VTK Vertical temperature eddy diffusivity (em“/sec) 


W Lead width (cm) 
jl Depth (cm) 
7A Surface ice thickness (cm3/cm@) 


A descriptive flow diagram of the program is shown in Figure 16. 
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Figure 16. Descriptive Flow Diagram of Program. 


Problem Initialization: 
DIMENSION/COMMON 
Function statements, 
Constants 


Input of Initial Conditions: 
Temperature TNEW, 
Salinity SNEN, 
Sigma-t SIGT, 
Equilibrium Ice Thickness ZI 


a 


Computation of Time-Dependent 
- Temperature and Salinity Fields 
TNEW, SNEW 





Computation of Surface Temperature Chances 
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Figure 16 (continued). 

















Computation of: 
Surface freezing temperatures TF, 

Surface temperature changes due to ice 
formation TIF, 

Iterative rate of ice formation HICE, 
Iterative volume of ice formed ZICE, 
Iterative salt rejected SSCHG, 

Resultant surface temperature and salinity 
fields (if melting occurred as a result of 
convective overturn, the resultant dilution 
is distributed over the mixed region of the 
water column) TNEW, SNEW 





Computation of Surface Ice Thicknesses ZI 
(iteration to account for freezing-point | 
changes aS a result of salinity changes) 


Vertical Stability Check a | | 
| OVRTRN | | 









Check for Surface Heat Brought 
From Depth Which Would Promote 
Melting 






Computation of: 
Rate of ice formation, 

Water-column heat-flux densities, 
Cumulative water column heat densities 





Fill Temperature and Salinity Arrays: 
I 
S 


Computation of Surface Heat-Flux Densities and 


Conversion to an Integrated Surface Temperature 
Gradient 
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